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结合 视觉 显著 模型 与 水 平 集 算 法 的 建筑 物 立 面 图 像 轮廓 快速 提取 


李 昌 华 ， 杜 文 强 ， 周 方 晓 


(西安 建筑 科技 大 学 信息 与 控制 工程 学 院 ,， 西安 710055) 


摘 要 : 针对 传统 轮 廊 提取 方法 对 建筑 物 立 面 图 像 处 理 存 在 的 处 理 耗 时 间 题 ， 提 出 一 种 结合 视觉 显著 模型 与 水 平 集 算 
法 的 建筑 物 立 面 轮廓 提取 算法 。 首 先 ， 通 过 视觉 显著 特征 与 超 像素 分 块 信 息 滤 除 大 量 无 关 背 景 ; 其 次 ， 对 所 得 的 显著 
分 块 ， 以 超 像素 分 块 为 单位 ， 据 建筑 物 的 纹理 特征 进行 算 选 ; 最 后 , 利用 水 平 集 算 法 对 所 得 初始 建筑 物 区 域 进行 演化 ， 
得 到 最 终 轮 廊 。 实 验 结果 表明 ， 较 传统 轮廓 提取 方 法 ， 在 处 理 效果 相当 的 情况 下 本 算法 处 理 速 度 得 到 有 效 提 升 。 
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Fast contour extraction from building facade image based on visual saliency model and level 
set algorithm 


Li Changhua, Du Wenqiangi, Zhou Fangxiao 
(School of Information & Control Engineering, Xi’An University of Architecture & Technology, Xi’an 710055, China) 


Abstract: Aiming at the time-consuming problem of traditional contour extraction method in building facade image processing, 
a building elevation contour extraction algorithm combined with visual saliency model and level set algorithm is proposed. First 
of all, using visual salient features and super pixel information to filter out the irrelevant background; secondly, it selects the 
significant blocks according to the building’s texture feature, with each Super pixel block as a processing unit; finally, using the 
level set algorithm optimizes the resulting building area to obtain the final contour. The experimental results show that the 
processing Speed of the proposed algorithm is effectively improved when the processing effect is equivalent, compared with the 
traditional contour extraction method. 
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和 直线 提取 的 方法 扩 匀 , 在 建筑 物 立 面 图 像 中 , 由 于 背景 复杂 ， 

及 成 像 过 程 中 引起 的 畸变 ， 各 种 边缘 检测 算法 无 法 准确 地 区 分 

图 像 型 三 维 重建 由 于 其 相对 低廉 的 成 本 ， 更 加 适合 于 大 规 。 全 部 建筑 物 边缘 ， 很 难 实现 全 部 轮廓 自动 提取 ; 9) 结合 主动 轮 

模 推 广 ,但 当前 图 像 型 三 维 重建 技术 依然 存在 处 理 效 率 不 够 高 、 ” 廓 模型 提取 目标 轮廓 的 方法 03-5。 主 动 轮廓 模型 9 (Snake 模 
精度 不 足 等 问题 。 若 能 在 三 维 重 建 的 前 期 ， 先 提取 其 外 轮廓 ， 型 ) 是 在 已 有 的 初始 轮 廊 曲线 的 基础 之 上 ， 赋 予 曲线 以 能 量 函 
针对 待 处 理 目 标 进行 预先 分 析 ， 将 这 些 信息 融合 到 现 有 的 三 维 ，” 数 ， 通 过 最 小 化 该 能 量 函数 ， 使 得 曲线 向 真实 轮廓 演化 ， 当 达 
重建 技术 中 ， 可 以 在 一 定 程度 上 提高 重建 的 效率 。 在 现实 生活 “到 能 量 最 小 时 ， 就 得 到 了 真实 轮廓 曲线 的 位 置 。Snake 模型 存 
中 ， 建 筑 物 轮廓 提取 和 对 提取 的 信息 进行 有 效 的 分 类 识别 在 在 实现 困难 ， 且 对 四 型 分 割 效 果 极 差 的 问题 ， 而 这 在 建筑 物 中 
市 化 规划 发 展 、 国 防 建设 、 抗 灾 减 灾 等 方面 有 着 重要 的 研究 经 常会 出 现 ， 所 以 Snake 模型 不 能 直接 使 用 。 而 主动 轮廓 模型 
义 ， 此 外 ， 对 于 图 像 检 索 、 三 维 模型 的 自动 生成 、 建 筑 物 识别 ”的 另 一 种 实现 方式 水 平 集 (LevelSet) 方法 上] 出现, 有 效 地 解决 
也 有 着 重要 指导 意义 。 其 中 建筑 物 轮廓 信息 提取 是 完成 三 维 建 ”了 此 问题 。 在 学 者 们 的 努力 下 ，LevelSet 方法 不 断 的 完善 ， 从 
筑 物 模 型 重 构 的 一 个 必要 条 件 ， 它 可 以 为 重 构 模 型 提供 更 多 的 Mumford-Shah 泛 函 ， 到 经 典 的 Chan-Vese 模型 ， 再 到 如 今 形 形 
场景 先 验 知识 和 信息 约束 条 件 。 目 前 常见 的 轮廓 提取 方法 主要 色色 的 LevelSet。 由 于 主动 轮 廊 模 型 主要 依赖 迭代 能 量 函 数 实 
有 三 类 : a) 预 先 处 理 一 批 常 见 的 包含 目标 的 模板 ， 利 用 先 验 目 岗 ， 其 处 理 效 率 由 初始 轮廓 和 人 烛 代 次 数 共 同 决定 ， 因 此 ， 一 般 
标 信息 ,采用 模板 匹配 的 方式 提取 目标 轮 廊 131, 此 法 实现 简单 ， 不 单独 使 用 。 
但 提取 建筑 物 有 诸多 限制 ， 且 提取 精度 有 限 ; b) 结 合 边 缘 检 测 针对 建筑 物 立 再 
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像 ， 上 述 方法 在 各 个 层面 各 具 优 和 劣 ， 但 
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普遍 存在 提取 效率 不 高 、 
图 像 分 块 的 基础 上 引入 视觉 显著 模型 ,提出 一 种 以 图 像 
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位 A 


两 类 : 


代表 


b) 针 对 特定 目标 ， 基 于 其 
模型 。 由 于 在 建筑 物 立 画 
E 以 及 显著 的 纹理 
以 建筑 物 区 域 为 视觉 搜索 目标 ， 以 区 


的 形 
显著 
和 纹 
并 将 
过 度 


纹理 


性 注意 


建筑 物 江 本 


经 实验 证 明 ， 该 方法 很 好 的 解决 了 建 多 


提取 问题 。 


耗 时 较 长 的 问题 。 为 解决 此 问题 ， 本 


以 超 像素 分 块 为 处 理 单位 ， 结 合 水 平 集 演化 算 
图像 轮廓 快速 提取 算法 ， 其 整体 框图 如 图 
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图 1 本 文 


法 整体 框图 


初始 轮廓 提取 
显著 图 的 获取 


人 类 视觉 注意 机 制 的 作用 在 于 从 复杂 的 视觉 场景 中 快速 定 


感 兴趣 目标 的 位 置 。 


为 Itti 模型 09]， 当 


a) 基 于 低级 视觉 特征 


当前 模拟 视觉 注意 的 计算 模型 主要 分 为 
的 自 底 向 上 的 视觉 计算 模型 ， 其 典型 
背景 分 布 较为 杂乱 时 处 理 结果 不 够 理想 ; 
高 级 视觉 特征 的 自 顶 向 下 的 视觉 ; 


图 像 中 ， 建 筑 物 区 域 具有 典型 


状 特 程 


特征 的 特性 ， 因 此 ， 本 文 利 用 视 沉 


模型 ， 
理 特 征 


作为 视觉 注 


其 作为 判断 该 区 域 是 
量 该 指标 来 确定 建筑 物 区 域 的 位 


或 的 形状 特征 
主意 模型 中 底层 特征 ， 计 算 区 域 显著 性 ， 
否 为 建筑 物 区 域 的 评价 指标 ， 最 终 通 


为 凸显 建筑 物 区 域 纹 理 


的 区 分 度 ， 文 献 [ 


并 提高 建筑 物 区 域 纹理 与 其 他 干扰 
20] 依 据 建 筑 物 区 域 纹 理 的 排列 规则 及 朝 


C 期 ] 
守 虽 华 ， 等 ， 结 合 视 党 显著 模 型 与 水 平 入 算法 的 过 各 扬 立 轴 让 衣 愉 订 记 
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玛 乘积 , YUo) 为 TC1o) 的 共 因 ，sgn 为 符号 
分 别 为 QDCT 变换 以 及 QDCT 反 变 换 。 
1.2 显著 图 像 分 块 
单个 像素 不 足以 表征 


函数 , QDCT 和 1QDCT 


建筑 物 立 面 的 特征 信息 ， 从 而 需要 对 
z 面 图 像 进行 分 块 处 理 。 为 加 快 后 续 处 理 过 程 ， 本 文采 
用 超 像素 分 块 与 视觉 显著 图 相 结合 的 方法 PC34， 在 SLIC 超 像 素 
分 块 算法 的 基础 上 ， 结 合 得 到 的 显著 图 ， 对 图 像 有 偏 分 块 ， 对 
显著 区 域 生 成 密集 超 像素 ， 在 非 显著 区 域 生 成 相对 稀 纹 的 超 像 
在 减少 全 局 超 像 素数 目的 同时 保留 了 边缘 细节 信息 ， 有 利 
于 后 续 的 提取 处 理 。 

通过 阔 值 化 显著 图 $_Map， 得 到 显著 区 域 SR 和 非 显 车 区 域 
NSR。 在 此 基础 上 ， 超 像素 分 块 算法 步骤 如 下 : 

a) 网 格 初始 化 。 对 显著 区 域 SR 和 非 显 著 区 域 YSR 分 别 设置 
间距 为 $1 、52 的 网 格 ， 初 始 超 像 素 中 心 选取 网 格 中心 ， 并 将 SR 
区 域 的 网 格 中 心 显著 性 设 为 0，NSR 区 域 网 格 中 心 显著 性 置 1。 


1 
5 时 (2 (3) 
$: 
5z= 芋 (人 


其 中 , N 表 示 图 像 总 像素 数 , c 表 示 显 著 像 素数 , t 表 示 非 显著 区 
域 NSR 与 显著 区 域 SR 超 像素 大 小 的 比值 ，K 为 超 像素 分 块 总 数 
目 。 

b) 一 步 聚 类 。 在 定义 的 距离 函数 Dij 下 ,对 每 个 超 像素 中 心 
Cx 对 应 的 搜索 区 域内 的 像素 进行 一 次 距离 计算 ， 并 将 该 像素 划 
分 到 最 近 的 超 像素 分 块 中 。 每 次 的 搜索 范围 为 以 其 所 在 超 像 素 
中 心 为 中 心 的 矩形 区 域 ， 搜 索 范 围 的 大 小 依据 超 像素 中 心 的 显 
著 性 而 不 同 , 若 中 心 显著 性 为 0, 搜索 范围 设 为 251*251; 否则 ， 
搜索 范围 设 为 25,*25,。 


dt 


向 一 致 的 特点 ， 结 合 非 下 采样 轮廓 波 变 换 (nonsubsampled Dix= (+ 人 mw) 
contourlet transform, NSCT) 的 多 尺度 、 多 方向 及 平移 不 变 等 特 (s_Map())? 
性 ， 将 NSCT 分 解 后 的 各 尺度 和 方向 子 带 的 局 部 纹理 能 量 作为 0 . 
其 纹理 特征 ， 再 结合 最 大 化 各 特征 的 信息 量 ， 得 到 图 像 的 显著 。 ”其 中 ，d。 表 示 像 素 ! 与 超 像 素 中 心 Ci. 之 间 的 颜色 距离 ,ds 表示 像 
区 域 。 但 由 于 该 文献 的 处 理 对 象 为 高 分 辩 率 遥感 影像 ， 相 比 于 。” 素 与 超 像 素 中 心 C 之 间 的 空间 距离 , S 表 示 网 格 大 小 ,S_Map(i) 
背景 更 加 复杂 的 建筑 物 立 面 图 像 ， 还 需 加 入 更 多 全 局 性 质 的 特 ”表示 图 像 中 像素 的 显著 性 。 
征 。 针 对 此 问题 ， 文 献 CI 以 显著 区 域 " 全 局 孤立 ， 局 部 紧凑 "为 虽 边 缘 细 化 。 遵循 SLIC 算法 框架 ， 查找 超 像素 中 心 Ci 搜索 
基础 ， 将 显著 目标 检测 问题 转换 为 马尔 科 夫 随机 游 走 问题 ， 提 ”范围 的 边缘 像素 P, 计算 像素 与 超 像素 中 心 Ci 之 间 的 距离 Dix， 
出 了 一 种 基于 混合 图 上 随机 游 走 的 显著 目标 检测 算法 。 但 由 于 当 Dix 大 于 一 定 的 阔 值 时 ,对 像素 户 重 新 分 配 最 近 的 超 像素 中 心 
其 构建 的 是 全 连通 的 混合 图 , 算法 时 间 复杂 度 为 om， 当 输 入 C，v。 随 着 迭代 过 程 的 进行 ， 超 像素 边缘 不 断 地 细 化 ， 直 到 超 像 
图 像 分 辩 率 较 大 ， 分 块 较 多 时 ， 并 不 适合 于 快速 提取 显著 区 域 。 素 目标 误差 函数 收敛 。 
的 情况 。 1.3 ”建筑 区 域 提取 

为 快速 计算 显著 特征 ， 本 文选 用 相位 谱 法 中 的 QDCT 对 于 获取 的 建筑 物 区 域 ， 一 方面 ， 其 中 可 能 参 杂 部 分 非 奸 


(quaternion DCT) 法 0 提取 图 像 的 显著 图 。 


其 中 : 


Sopcr (lo) 9g* [T(10) a T(10)] 
T(10) = IQDCT!(sgn(QDCT*(10))) 


(D) 
2 


5 表示 提取 的 显 


著 图 ，g 表 示 高 斯 滤波 器 ，lo 是 在 原始 图 


像 1 各 通道 图 的 基础 上 


转变 而 来 的 图 像 矩 阵 ， 运 算 符 。 表 示 哈 达 


筑 物 区 块 ， 如 树木 、 草 地 等 ， 另 一 方面 ， 该 区 域 可 能 并 非 包含 
完整 的 建筑 物 区 域 。 因 此 ， 需 要 对 显著 分 块 采取 进一步 筛选 处 
理 。 

建筑 物 立 面 图 像 具 有 “形状 规则 , 纹理 明显 ”的 特点 。 其 中 ， 


“形状 规则 ? 指 典型 的 建筑 物 立 面具 有 规则 的 几何 形状 棱角 ， 
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常 以 直角 和 矩形 为 主 ， 几 何 形态 多 呈现 为 规则 的 多 边 形 结构 ， 而 
“纹理 明显 ” 指 建筑 物 立 面 含 有 明显 的 纹理 特征 ， 且 纹理 分 布 通 
定 的 方向 性 ， 多 体现 为 与 建筑 物 轮廓 方向 近似 一 致 或 
而 超 像 素 分 块 中 可 以 很 好 的 六 理 特性 ， 因 此 ， 
本 文选 取 纹 理 特征 作为 进一步 筛选 图 
现 有 的 纹理 特征 提取 方法 ， 
模型 法 、 信 号 处 理 法 以 及 结构 法 R31， 而 其 中 统计 法 由 于 其 思 
容易 实现 而 广泛 使 用 ， 本 文选 用 统计 法 中 的 灰 度 共生 条 
阵 (GLCM) 法 提取 建筑 物 纹 理 特征 。 为 了 更 好 的 描述 分 块 特 
性 ， 本 文 在 特征 参数 中 加 入 分 块 的 颜色 分 量 ， 因 此 ， 提 取 分 块 
的 特征 向 量 F = (有 i, 所 , …, 户 ) 共 有 7 个 参数 。 
1) 颜色 分 量 (color components) 

fe= Bi Lli,)) (7) 
其 中 ，N 表 示 分 块 内 总 灰 度 级 别 数 ，c = 二 2, 3 表示 不 同 的 颜色 
通道 ，I 表 示 分 块 图 像 1 的 通道 ，i，j 表 示 不 同 的 灰 度 级 别 。 
2) 纹理 


Er 


十 止 父 。 


以 分 为 四 类 : 统计 法 、 
思想 


主 安 品 


简单 ， 


有 反差 (texture contrast) 
hh =i —)2P(i,)) 
k 生 和 矩阵 基础 上 提取 的 


(8) 
k 生 概率 矩阵 。 纹 


其 中 , P 为 在 分 块 灰 度 


办 
像 超 像素 分 块 的 主要 依据 。 


ChinaXiv 合 作 期 刊 


李 昌 华 ， 等 : 结合 视觉 显著 模型 与 水 平 集 算法 的 建筑 物 立 面 图 像 轮廓 快速 提取 


fe = 


—ZijP(i,)) :logz Pli,)) 


(10) 


纹理 箭 值 表明 了 纹理 灰 度 分 布 的 复杂 程度 ， 箭 值 越 大 ， 纹 


健 
Lal 
读 


里 相关 性 (texture correlation) 


f= zm LY PD) — pry) 


其 中 ，j，by 分 别 表示 菜 一 分 块 的 共生 概率 
直方 向 的 均值 , ox 
纹理 相关 性 在 一 定 程度 上 反映 了 分 块 图 像 中 


灰 度 分 布 越 均匀 ， 相 关 性 值 越 大 。 


ay 分 别 表示 水 平方 向 和 竖 直 方向 的 标准 


4D) 
矩阵 水 平方 向 和 坚 
差 。 
像素 间 的 灰 度 关系 ， 


首先 ， 建 筑 物 立 面 图 像 中 背景 主 
显著 分 块 以 各 颜色 分 量 均值 为 特 得 
类 的 聚 类 中 心 分 块 Ct，C2:， 随后 ， 
著 分 块 进行 盘 选 ， 剔 除 其 中 属于 背景 分 块 
筑 物 立 面 图 中 干扰 信息 多 有 分 布 集中 的 特点 
则 多 体现 为 以 显著 区 域 边界 开始 的 某 几 个 连 
过 搜索 显著 分 块 边界 中 与 C1 或 6 颜色 特征 相 
这 些 分 块 为 中 心 搜寻 其 邻 域 中 与 其 相似 的 分 


E 进 行 一 


岂 的 基础 上 ， 以 显著 区 域 边界 的 超 像素 分 块 
或 分 块 搜索 ， 以 分 块 特征 向 量 F 之 间 的 欧式 
， 往 选 建筑 分 块 ， 排 除非 建 


量 标准 筑 元 素 ; 


区 域 即 可 得 到 建筑 物 初始 轮廓 。 


初始 轮廓 


了 
是 


理 反 差 反 映 图 像 的 清晰 度 和 纹理 沟 纹 深浅 的 关系 ,反差 值 越 大 ， 
纹理 越 清 晰 。 
3) 纹理 能 量 (texture energy) 
f= Pi,)) (9) 
纹理 能 量 反 映 图 像 灰 度 分 布 均匀 程度 和 纹理 粗糙 细 读 ， 大 
的 能 量 值 表征 一 种 较 均匀 和 规则 变化 的 纹理 模式 。 
显著 
分 块 
非 显 著 两 类 、| 中 心 分 块 
分 块 聚 类 CuC2 


图 2 建筑 区 域 提取 框图 


2 ”水 平 集 轮廓 优化 


] 水 平 集 方法 的 图 像 轮廓 优 化 主要 是 将 平面 闭合 曲线 表 
将 二 维 空间 的 曲线 演化 问题 转变 
为 三 维 空间 的 曲面 演化 问题 ， 通 过 最 小 化 定义 的 能 量 函 数 获取 
最 优 轮廓 。 本 文 用 于 建筑 立 面 轮廓 优化 的 算法 是 在 水 平 集 算法 
的 一 种 实现 方式 DRLSEP“55 的 基础 上 稍 作 改动 实现 的 .DRLSE 
水 平 集 算法 无 须 重 新 初始 化 水 平 集 函 数 ， 在 很 大 程度 上 提高 了 


全 


采 
示 为 水 平 集 函 数 的 零 水 平 集 ， 


P(9) =ih (v8|— D2dxdy 


央视 为 背景 分 块 ， 从 显著 分 块 中 滤 除 ， 最 后 ， 


因此 将 
得 到 


次 两 类 聚 类 ， 


以 分 块 C1 和 C, 为 参照 ， 对 


的 部 分 ， 考 虑 到 建 
， 而 在 显著 分 块 中 
续 区 域 ， 因 此 ， 通 
似 的 分 块 ， 然 后 以 
块 ， 将 搜索 到 的 分 
在 得 到 的 显著 分 
为 起 始 ， 对 其 8 邻 
距离 作为 相似 性 度 
最 终 ， 合 并 建筑 物 


(13) 


外 部 驱动 能 量 Eg(9$) 的 定义 如 式 (14)~(16) 所 示 , 其 中 ,Lg(9$) 
表示 零 水 平 集 的 加 权 弧 长 , Ag(8) 表 示 零 水 平 集 的 加 权 面 积 , 参 


数 4 和 a 为 对 应 惩罚 系数 。 


Eg(@) = ALg(0) + Ag(9) 
Lo(@) = fo glo)6e(9) VOlaxdy 
Ag($)= fo g(a)He(—$)dxdy 


44) 
(15) 
(16) 


水 平 集 演化 效率 , 其 能 量 E($) 定 义 为 内 部 约束 能 量 PC9) 与 外 部 ”其 中 ，65 为 标准 差 为 0 的 高 斯 函数 ，H。， 5 分 别 为 Heaviside 函 
驱动 能 量 妨 (的 两 部 分 ， 如 式 (12) 所 示 。 数 和 Dirac 函数 , g(lo) 表 示 图 像 边缘 指示 函数 , 其 定义 如 式 (17) 
E($) = MP(9) + Eg(®) (12) ”所 示 。 

其 中 ，$ 二 《gry} 表示 水 平 集 函 数 ， P(g) 保 证 水 平 集 函数 在 i 本 

演化 过 程 中 保持 为 一 个 符号 距离 函数 ，/ 为 对 应 大 于 零 的 惩罚 到 

系数 ，B (内 驱 动 水 平 集 朝 着 目标 边缘 运动 。 由 于 建筑 物 立 面 图 像 中 存在 许多 纹理 细节 ， 水 平 集 在 演化 
内 部 约束 能 量 P( 风 定义 如 式 (13) 所 示 ， 是 对 水 平 集 函 数 偏 。 过 程 中 容易 受 其 影响 ， 同 时 ， 建 筑 物 立 面 真实 轮廓 多 处 于 超 像 

离 符 号 函数 的 度量 ，0 表 示 整 个 图 像 域 ，V 为 梯度 算 子 。 素 分 块 边界 处 ,， 因此， 本 文 在 DRLSE 外 部 驱动 能 量 本 (内 的 


合作 期 刊 | 


Chinaxn 
录用 稿 李 昌 华 ， 等 : 结合 视觉 显著 模型 与 水 平 集 算 ; eA 决 如 提取 


础 上 ， 引 入 分 块 能 量 项 Cy(8)， 如 式 (18)(19) 所 示 ，Co ($8) 用 于 加 立 面 图 像 的 部 分 边界 丢失 , 本 文选 取 2.1 的 折 中 方案 ; 参数 y 的 
快 处 理 超 像素 分 块 ， 在 超 像 素 分 块 中 增强 水 平 集 函数 值 ， 而 在 。” ” 取 值 不 宜 过 大 ， 且 应 小 于 参数 a， 本文 取 y=0.3。 水 平 集 各 参数 
超 像 素 边界 处 削弱 水 平 集 函数 值 ， 保 证 最 终 的 目标 轮廓 尽 可 能 ” 取 值 情况 如 表 1 中 所 示 。 


不 受 建筑 物 立 面 细节 影响 。 表 1 水 平 集 轮廓 优化 主要 参数 设置 表 
Eg(@) =415(0) +a4o(b) 二 7CrC0) (18) 参数 At 卫 A a y 


Cg($) = fo gM He(—p) dxdy (19) 取 值 5 02At 5 21 03 
其 中 : M 为 与 图 像 1 同 维度 的 超 像 素 分 割 掩 码 和 矩阵，g(L) 表 示 图 


像 的 超 像素 边缘 指示 函数 , 其 定义 如 式 (20), g(L) 在 超 像素 边缘 图 像 127788 的 水 平 集 轮廓 优化 实验 结果 如 图 4 中 所 示 ， 
处 取得 较 小 值 。 其 中 ， 图 4-a 为 初始 轮廓 ， 图 4-b 为 30 次 迭代 后 的 轮廓 ， 图 4- 
RE 而 c 为 50 次 迭代 后 的 轮廓 。 可 以 看 出 ， 随 着 欠 代 次 数 的 增加 ， 图 
elven 像 轮 廓 逐渐 趋 于 其 真实 轮廓 ， 对 于 底部 不 太 明 确 的 边界 ， 也 达 
利用 梯度 下 降 法 求解 式 (18) 的 极 小 值 , 等 价 于 求解 式 (2D) 的 “到 很 好 的 收敛 。 
稳 态 解 


= et -div (WE)) + Ga (90) BE) + 


aos(b)g(Wc) +y6(P)g(Lo) (21) 
将 式 (21) 离 散 化 ， 并 移 项 可 得 水 平 集 能 量 演化 式 ， 如 式 (22) 
所 示 。 


(a) 初 始 轮廓 (b)30 次 迭代 (0)50 次 锡 代 


人 + AtL (22) 
= bey (Px) 图 4 水 平 集 轮廓 优化 结果 (At=5，o=2.1， 人 0.3) 


其 中 , At 为 大 于 零 常数 , 表示 演化 步 长 , L(9#y) 为 式 (21) 的 右 端 


己 
= 
o 


通过 前 述 的 处 理 步 又 ,已 获取 到 了 有 关 建 筑 物 的 大 致 轮廓 ， 3 ”实验 与 分 析 

因而 水 平 集 方法 处 理 不 需要 做 全 局 的 处 理 ， 选 用 罕 带 水 平 集 方 本 文 的 实验 环境 为 iMac, Intel Corei5 处 理 器 , 主 频 3.4GHz， 
法 即 可 实现 ， 水 平 集 轮廓 优化 流程 如 图 3 中 所 示 。 8GB 内 存 ，NVIDIA GeForce GTX 775M 显示 卡 ，2GB 显存 ， 
软件 环境 为 MATLAB R2016a。 实 验 数据 主要 取 自 MSRA10K 


数据 集 。 图 5 为 本 文 算法 的 轮廓 提取 效果 展示 图 ， 其 中 由 左 向 

依据 水 平 集 函 数 堆 右 依次 为 原始 图 像 (图 $ (a))、 显 著 分 块 图 (图 $ (b))、 经 过 

交叉 更 新 窄带 B 建筑 区 域 提取 的 分 块 (图 5 (c))、 水 平 集 优化 后 的 最 终 轮 廊 图 

(图 5(d)) 以 及 采用 gPb-ucm 法 P9 提 取 的 轮廓 图 (图 $(Ce))。 
平 集 函 数 中 罕 带 

”8 等 参数 初始 化 图 中 蓝 绿色 的 线条 表示 超 像 素 分 块 边界 ， 而 红色 线条 表示 最 终 

的 轮廓 ， 可 以 看 出 本 文 的 轮廓 提取 算法 提取 的 轮廓 基本 收敛 于 


窄带 B 连 续 无 变化 
或 迭代 次 数 达 上 限 


建筑 物 的 真实 轮廓 处 。 

本 文采 用 精确 率 (precision)、 召 回 率 (recall) 和 FF 值 (F- 
measure) 评 价 提 取 轮 廊 的 质量 。 精 确 率 P 表征 检测 出 的 信息 正 
确 部 分 的 比例 ， 召 回 率 R 表示 有 效 信息 被 正确 检索 的 比例 ， 其 
计算 公式 如 式 (23)(24) 所 示 。 其 中 ，TP 表示 正 类 被 判断 为 正 类 
的 情况 ，FP 表示 负 类 被 判断 为 正 类 的 情况 ，FN 表示 负 类 被 判 
断 为 正 类 的 情况 。 精确 率 和 召回 率 的 取 值 位 于 0-1, 取 值 越 接 近 

水 平 集 轮廓 优化 的 主要 参数 有 At，h, 入 ，a 和 Yy。 参 数 At 用 ”1,， 表示 算法 性 能 越 好 , 但 二 者 通常 存在 矛盾 的 情况 ， 此 时 ， 通 
于 控制 曲线 演化 速度 ， 由 于 是 在 前 期 得 到 初始 轮廓 的 情况 下 进 ” 过 引入 三 者 的 加 权 平 均值 F 值 进行 评价 。F 值 的 计算 公式 如 式 
行 演化 , 因此， 可 以 将 At 设置 大 一 些 ， 本文 取 值 为 5; 参数 为 ” (25) 所 示 ， 其 中 a 为 引入 参数 。 当 a 取 值 为 1 时 计算 的 F 值 即 为 
符号 距离 调整 权重 ， 文献 PC9 中 指出 4 与 At 的 乘积 需 满足 4* Fl-measure。 


依据 窄带 B 更 新 水 平 集 


芭 带 水 平 集 轮 廓 优化 流程 


图 3 


rr] 


At < 1/4， 本 文 几 取 值 为 02/At; 参数 入 ，a, y 分 别 为 外 部 驱动 P= (23) 
能 量 项 各 分 量 权重 ， 由 于 DRLSE 对 参数 上 和 》 不 敏感 ， 本 文 R= (24) 
》 的 取 值 沿用 文献 Pg 中 取 值 = 5， 非 零 的 参数 a 会 导致 优化 的 ho, — Catern 三 
轮廓 与 真实 轮廓 发 生 一 定 的 偏差 ， 而 过 大 的 a 又 会 导致 建筑 物 th) 
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5 本文 算法 轮廓 提取 效果 


从 表 2 中 可 以 看 出 ， 本 文 算法 与 文献 [26] 的 轮廓 提取 方法 
相 比 ， 精 确 率 稍 有 提升 ， 而 召回 率 略 有 降低 ， 在 Fl 值 整体 相 
当 的 情况 下 ， 耗 时 缩减 了 近 一 半 ; 相 比 于 grabcut 算法 ， 由 于 
引入 超 像素 和 纹理 信息 ， 精 确 率 和 召回 率 都 略 有 提升 。 


表 2 算法 性 能 参数 对 比 


| 


| 


方法 耗 时 /s ”精确 率 ” 召回 率 Fl 
gPb— 

n 4.77 82. 7% 92. 5% 87. 3% 
Ucm' 
grabcut 3. 29 79. 2% 86. 1% 82. 5% 
本 文 方法 2. 56 85. 9% 89. 6% 87. 1% 


本 文 结 合 视觉 显著 图 与 图 像 分 块 思想 ， 利 用 视觉 显著 图 吻 
除 无 关 背 景 ， 在 显著 图 分 块 的 基础 上 ， 结 合 建筑 物 特有 的 颜 
色 、 纹 理 以 及 形状 特征 ， 过 滤 树 木 等 干扰 ， 进 而 融合 建筑 物 区 
块 ， 得 建筑 物 初始 轮廓 ， 最 后 结合 水 平 集 算法 演化 ， 得 到 经 过 
优化 的 建筑 物 轮廓 。 通 常 算法 处 理 精 度 与 提取 速度 之 间 存 在 矛 
盾 性 ， 高 的 精度 ， 意 味 着 相对 较 低 的 处 理 速度 。 实 验 表 明 ， 本 
文 的 轮廓 提取 算法 对 于 一 般 的 建筑 物 立 面 图 像 可 以 快速 地 提取 
其 轮廓 ， 并 且 可 以 通过 调整 水 平 集 算法 的 迭代 次 数 进而 保持 轮 
廓 精度 与 提取 速度 之 间 的 平衡 。 同 时 本 文 的 建筑 物 立 面 轮廓 提 
取 算 法 也 存在 不 足 之 处 ， 主 要 针对 单 目标 情况 进行 处 理 ， 而 算 
法 得 处 理 效率 也 依赖 初期 显著 图 ， 当 显著 区 域 偏 离 真实 建筑 物 
区 域 时 ， 处 理 效率 不 够 理想 。 
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